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1. Introduction: 

Relativistic aspects of quantum entanglement, especially in Einstein, Podolsky and 
Rosen (EPR) correlations [T] and the Bell inequality [21 [31 HJ 13 El [7J, have recently 
attracted much attention. The quantum states in inertial frame of reference transform 
under Wigner rotation [8] via Non-Abelian continuous group SU(2) representation in 
momentum space. Therefore, under Lorentz transformation, spin becomes entangled 
with the particle's momentum. Photon polarization qubits behave similarly and linear 
polarization states of photon will be seen in moving frame in the form of depolarized 
states [9]. 

Other aspects of relativistic entanglement are based on spinor formalism, which is 
designed for four-dimensional space-time [10]. In Refs. [TT| IT2] \13\ [TH [15] . Czachor 
and co-worker have used spinor method and shown that the helicity-momentum states, 
which are the irreducible representation of Poincare group, are the projection of the 
relativistic spin operator on timelike direction in momentum representation. 

Separability of the quantum states and quantification of entanglement in composite 
systems are perhaps the most important features of quantum information. The first 
related criterion for distinguishing entangled states from separable ones, is the positive 
partial transpose (PPT) criterion, introduced by Peres [16] . Nonetheless, the strongest 
manner to characterize entanglement is using entanglement witnesses EWs [T71 [IB] . 
An EW is an observable W whose expectation value is nonnegative on any separable 
state, but strictly negative on an entangled state. Recently, there has been an increased 
interest in the NLEWs because of their improved detection with respect to linear EWs. 
An NLEW is any bound on nonlinear function of observables which is satisfied by 
separable states but violated by some entangled states [HI [20]. For the first time 
in Refs. [22, [23] , it was actually shown that EW can be expressed as a measure of 
entanglement. However, the main objective of this paper is to show that NLEWs can 
also be very helpful to the quantification of entanglement. 

In our previous work [24] , we studied spin- momentum correlation in single-particle 
spin- 1 quantum states by using concurrence and showed that the amount of spin- 
momentum correlation depends on the angle between spin and momentum. But here, 
We quantify the spin-momentum correlation for a massive spin half and one relativistic 
single-particle in 2 £g> 2 (as two-qubit system) and 2 <g) 3 dimension Hilbert space using 
the NLEWs. For simplicity, instead of superposition of momenta, we have used two 
momentum eigenstates (pi and p 2 ) and in 2D momentum subspace, our results suggest 
the effect of the Lorentz transformation would decrease both the amount and the region 
of entanglement. 

This paper is organized as follows: Sec. II, is devoted to single-particle spin half and one 
quantum states. In Sec. Ill the spin-momentum correlation is calculated in relativistic 
single-particle spin half and one mixed states by using the NLEW. Finally, in Sec. IV 
the results are summarized and conclusions are presented. 
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2. Relativistic one massive spin-| and one particle quantum states in 
helicity basis 

2.1. Unitary representation of Poincare group 

The Poincare group ISO(3,l) is a semidirect product of S0(3,l) group of Lorentz 
transformation with the following generators 

P»,Ji" fjL,u = 0,1, 2, 3 (1) 

where P^,fi = 1,2,3 are components of momentum operators, and P° = p° = 
y/rnFc 2 + p 2 is Hamiltonian, and J = {J 23 , J 31 , J 12 } are components of angular 
momentum and K = {J 10 , J 20 , J 30 } are boost operators. Likewise, unitary 
representation of Poincare group can be written in the following form: 

K = -^V p --P^, (2a) 
mc + p u 

J = -ip x V p + s, (26) 

P = p, (2c) 

where Vp is the partial differential as ^ and s denotes the finite dimensional angular 
momentum corresponding to s 2 = s(s + 1). We know that the Poincare group is of rank 
2, so there are only two independent Casimir invariant operators, which are squared 
mass and the square of the Pauli-Lubanski vector W M that commute with all generators 
of the algebra. Let's use this definition to write down the form of Casimir's operators. 

the Pauli-Lubanski vector can be written as: 

W„ = \e^ P aJ vp P a (4) 

where e^ upa is anti symmetry Levi-Civita tensor and takes the value +1 if /jupa is an 
even permutation, -1 if it is an odd permutation and zero otherwise ( e i23 = — e 0123 = 1) 
. It is common to denote the Pauli-Lubanski vector as 

W p = {W°, W) = (P ■ J, P J — P x K). (5) 

We can define the popular center-of-mass position operator as 

Q = ^K + K-I]. (6) 

^ -M) -M) 

Then, the orbital angular momentum is 

L = Q x P. (7) 
Therefore, We can write a relativistic spin operator in the following way 

S = J - L, (8) 
and using the Pauli-Lubanski vector, we get 

o W /n\ 
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So, the spacelike component of Pauli-Lubanski vector is proportional to the relativistic 
spin operator. If we choose the timelike vector b a = (b°, 0, 0, 0), then after projection of 
Pauli-Lubanski vector on it, we obtain 

b a W a = b°W° = b°p ■ J, (10) 

which is helicity. On the other hand, the projection of Eq.Q onto the vector b a = (0, b) 
have eigenvalues as 

±1 



\± = —^(p-by + m^. (11) 
If the world- vector b a could be a null tetrad, i.e., b 2 = b a b a = 0, we get 

A± = ±i(p-6). (12) 

In momentum representation, we can write two energy-momentum world- vectors as 
a linear combination of the null tetrad and after projection of RSO onto the null tetrad 
(see Eq. ffT2]) ). we get the following eigenvalue which is similar to eigenvalue Pauli's 
matrix, i.e., a z ( for details see [H] ). 

A± = ±\. (13) 
2.2. single-particle spin-^ quantum state 

Here in this paper, quantum state is made up of a single-particle having two types of 
degrees of freedom : momentum p and spin a. The former is a continuous variable with 
Hilbert space of infinite dimension while the latter is a discrete one with Hilbert space 
of finite dimension. For simplicity, instead of using the superposition of momenta, we 
use only two momentum eigenstates (p% and p?). However, we restrict ourselves to 2D 
momentum subspace with two eigenstates pi and p 2 , so the pure quantum state of such 
a system can always be written as 

2 

=E E CijlPua), (14) 

1=1 (7=0,1 

where |pi(2)) and \a) are eigenstates of momentum and spin operators, respectively, c^- 
are complex coefficients such that Yhj \ c ij\ 2 = 1- 

A bipartite quantum mixed state is defined as a convex combination of bipartite pure 
states ([HI) , i.e., 

Pi = EWi>M, (15) 
i=i 

where Pj > 0, J2i Pi — 1- The subscript " | " refers to the spin-| and |^) (i = 1,2,3,4) 
as four orthogonal maximal entangled Bell states (BD) belong to the product space 
% p ® %cr and are well-known as 

|^i> = ^(bi,0) + b 2 ,l>) ; (16a) 
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|V>2) = ^(bi,o)-b2,i>), (166) 
|^> = ^(ba,0) + b 1 ,l», (16c) 

|^4> = -^(|pa,0>-|p 1 ,l». (16eO 

in terms of momentum and spin states. The kets |0) and |1) are the eigenvectors of spin 
operator a z . We assumed that spin and momentum are parallel in the z-direction and 
in this case, the single-particle spin-^ state can be considered as a two-qubit system. 
For an observer in another reference frame S', is described by an arbitrary boost A in 
the x-direction. The transformed BD states are given by (see Appendix A), 



hk) — > u(AM) = J^Y, D MW(A, P )M), (i7) 

VP c t 

where U (A) is a unitary representation of Lorentz transformation. It can be calculated 
that will be orthogonal after Lorentz transformation, i. e., 

(A^|A^>=%. (18) 

The BD density matrix ffl5l) . which describes the state of the single-particle at non- 
relativistic frame, is exchanged to the density matrix pf after Lorentz transformation, 

2 

i.e., 

pi — ► U(A)p h U(A)\ (19) 

therefore 

pi = U(A)piU(A)1 = Y,Pi\Wi)(W>i\- (20) 

2 2 -i 

i=i 

2.3. single-particle spin-1 quantum state 

We assumed that for spin one, the particle moves with two momentum eigen-state p\ 
and P2 along the y-axis. For an observer in another reference frame S' described by an 
arbitrary boost A given by the velocity v' in the z-direction, we have 



Pl(2) = (0,±p 1(2) ,0,£) = (0,m^V 1(a) 2 -l, (W 1(a) ), (21) 



where 



= ™ 2 .7' I(2) = f-^- ,y = jr^e- - c ( 22 ) 

According to Ref.[25], the Wigner representation for spin j is calculated as follows: 



V j , (W(A,p)) = e^ {5 -'^+^ u ' a){j+a+1)+5 -'---^ {j+CT){j ' CT+1) \ 



(23) 
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where 9 is angle around the x-axis. After some mathematical manipulations for spin 

one we get 



I cos0x(2)+l isin0 1(2 ) cos0 1(2 )-l \ 
' 2 72 2 ' 

1J21 COS 01(2) m 



72 """^J V2 

V 2 71 2 



(24) 



for simplicity, assume that p\ = —p2 = p, then 9i = — 6> 2 = 9 . 

In order to consider single spin one particle mixed state under Lorentz transformation 
we define the following density matrix in the rest frame as follows: 



1 — x — y 

Pi = ^ 1^5) (V's I +2/|^6)(^e| H 7. I,x + V < 1, 

o 



(25) 



The subscript "1" refers to the spin one and 1^5(6)) are maximally entangled pure states 
given by 

1 



|V%> = -^(|p 1 ,-l> + |pa,l», 



V2 



(bi,-i)-|P2,i)), 



(26 a) 
(266) 



where \pi(2)) are two momentum eigen states of particle and {| ± 1), |0)} are the three 
component spin ones as 



(27) 





(1\ 










ll> = 





,|o> = 


1 


,1-1) = 















K 1 ) 



By using (|24l) . we can obtain the relativistic density matrix (|25|) as 
p^p^ = U{K) Pl U\K). 



(28) 



3. Measure of entanglement of single-particle states using NLEW 

3.1. Entanglement Witnesses 

An entanglement witness W is an observable that reveals the entanglement of some 
entangled state p, i.e., W is such that Tr(W • p s ) > for all separable p s , but 
Tr{W ■ pen) < 0. The existence of an EW for any entangled state is a direct consequence 
of Hahn-Banach theorem [26] and the space of separable density operators is convex and 
closed. Geometrically, EWs can be viewed as hyper planes that separate some entangled 
states from the set of separable states and, hyper plane indicated as a line corresponds 
to the state with Tr[W ■ p s ] = 0. According to Refs. [27] 128] . an EW will be optimal if, 
for all positive operators P and e > 0, the operator 

W = (1 + e )W - eP, (29) 
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is not an EW. 

When talking about EWs, one has to take an important point into consideration: the 
so-called decomposable EWs (DEW), which can be written as 

W = P + QT A , V,Qi>0, (30) 

where the operator Qi is positive semidefmite. It can be easily verified that such 
witnesses cannot detect any bound entangled states. W is non-decomposable EW if 
it can not be put in the form (1301) (for more details see [29]). One should notice that 
only non-decomposable EWs can detect PPT entangled states. 

3.2. Measure of entanglement of single spin-^ particle using NLEW 

According to Ref. [30], [15] we present an NLEW for a bipartite system "H 2 ® % as 
follows 

4 

w 202 = h ® h + E AA ® 4h (si) 

where J 2 is a 2 x 2 identity matrix, A^j are some parameters, and OiS are Hermitian 
operators from first ( second ) party Hilbert space as following 

Oi = |0)(0|, (32a) 

O a = |l)(l|, (326) 

3 = ^=(0 1 + 2 ), (32c) 

04=^(01-02). (32 d) 

We introduce the maps for any separable state p s which map the convex set of separable 
states to a bounded convex region named as feasible region (FR). Then, recalling 
the definition of an EW, we imposed the first condition which is the problem of the 
minimization of expectation values of witness operators W with respect to separable 
states, i.e., 

min Tr{W ■ p.) > 0, (33) 

In the second step, for a given p ent , we imposed the second condition for an EW, 
Tr{\V ■ p en t) < 0. Now the objective function ( which will be minimized ) is Tr(W ■ p ent ), 
and the inequality constraints come from the first step solution. So, this problem can 
be written as a convex optimization problem. As mentioned above, we will use two 
steps towards finding the parameters Aij for the density matrix and fully characterize 
NLEWs based on exact convex optimization method for single spin-| quantum mixed 
state. Now we define pi as follows: 

2 

p~i . =Tr {pi -Oi^Oj), (34) 
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Obviously, these components are associated with components of witness matrix, i.e., 
A^ v . If matrix Z is a symmetric matrix (i.e. Z l = Z) then pi = —2AZ or 

Z= 1 -A t p h = 1 -(pVp l2 )K (35) 

where 

A=-\p\Z~\ (36) 
So, after some calculation we get 

An = A 22 = ~( - ' ^ ' ^ _ i), (37,, 







-P 2 


+ P3 + P4 






f Pi 


-Pa - P4) 2 


Ir- 


Pi4 


-p 2 . 


- Pa - P 4 



^12 = A11 = »( 7==^ = = - 1), (376) 
2 y(Pi + P 2 - Pa - P 4 ) 2 

, "Pi + P2 — P3 + Pi \ / Q „ s 

^33 = / ), (37c) 





-p 2 


+P3-P4) 2 


Pi- 


-p 2 - 


- p 3 + p 4 




-p 2 


-P3+P4) 2 



A 4A = 1 ' ' , (37d) 

^ - P 2 - P 3 + P 4 ) 2 

and other parameters are zero. Then, by using (I49p the NLEW can be written as 



>V 202 = h ® i" 2 + An(Oi ® Oi + 2 ® 2 ) + A 12 (d <8) £> 2 + C 2 g> Ci) 

+A33O3 ® £> 3 + A44O4 ® 4 - (38) 

Then measuring the observable W 282 of pi gives a good estimate of the expected value 
of W 202 , 

Tr[W 202 • pi] = ~{1 - |Pi + P 2 - P 3 - P 4 | - |Pi - P 2 + P 3 - P 4 | - |Pi - P 2 - P 3 + P 4 |}, 
2 z 

using X)i P?. = 1 we obtain 

Tr\Wwpi] = (l-2Pi). (39) 

We use the concurrence to measure the entanglement between spin and momentum of 
single particle. It is defined as 

C(p) = max(0, v\ — v<i — u 3 — z/ 4 ), 

where the quantities v\ > u 2 > V3 > v± in decreasing order are the square roots of 
the eigenvalues of the matrix 

R = p(<7 y <g> (y y )p*{o y <8> (Ty), 
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for Bell- diagonal density matrix of ( 11 5 p after some calculation we get 

C(pi) = (2P 4 - 1), (40) 

and Eq. (l39|) . we obtain 

C{pA = -Tr[W 2m -p h l (41) 



where Tr[W 2( g, 2 -pi] is zero when p = | and exactly coincides with the concurrence 
of two-qubit Bell-diagonal mixed state ( 1T51) . 

Moreover, we have made the NLEW for relativistic density matrix (120]) and taking the 
trace we obtain: 

Tr[W 202 ■ pi] = A — \(\/B-AVC + \fB + 4VC), (42) 

2 4 

where 

1 



A — 1 — -{sJ{Px + P2-P3- Pa) 2 + \]{Px -P2-P3 + Pa) 2 } cosfi p , (43a) 

B = 8(PiP 3 + P2P4) cos 2 n p - 2(P 2 + P 2 + P 3 2 + P 2 ) (-3 + cos 2Vt p ) , (436) 
C = 4(P + P 3 ) 2 (P 2 + P 4 ) 2 + 4((P 2 - P 2 ) 2 + (P| + P 2 ) 2 ) sin 2 fi p + 

4(P - P 3 ) 2 (P 2 - P 4 ) 2 sin 4 2n p . (43c) 
We suppose that the Wigner angle Q p is defined by (see Appendix A) 

^P2 (44) 
After some calculations and using ( )39l) . we can obtain 

Tr[W 202 • pi] < Tr[W 202 • pi] => C(p h (Q p = 0)) > C(pi(Q p )). (45) 

This result indicates that both the amount and region of entanglement decrease(see 
Figure. 2). 



3.3. Measure of entanglement of single-particle spin-one mixed state using NLEW 

We know that positive partial transpose (PPT) criterion is a necessary and sufficient 
condition for determining entangled states living in Hilbert spaces H 2 ®H 2 and H 2 ®H 3 . 
But, with using the NLEW we want to measure of entanglement of bipartite system in 
H 2 <S> P3 Hilbert space. 

In this section, we want to calculate the geometric measure of entanglement of 
single spin one relativistic particle mixed states in the 2 <g> 3 Hilbert space which are 
spanned by the following basis vectors 

{bi( 2 )>® |±l»i( 2 )>® |0>}, (46) 
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We consider the density matrix of (I25p and obtain the FR by imposing PPT conditions 
with respect to each party. The PPT condition was applied on the first party and the 
eigenvalues of p T in the rest frame (i.e., = 0) are given by 

\ 1 = hl + 2x-4y), (47a) 
o 

\ 2 =l(l-x-y), (475) 
o 

\ 3 = ±(l-4x + 2y), (47c) 

\ 4 = hl + 2x + 2y), (47 d) 

o 

The positivity of density matrix ( 1251) imposes the following constraints on the parameters 
\ 1 = hl-x-y), (48a) 



X 2 = hl + 5 x-y), (48b) 


A 3 = Ul-x + 5y), (48c) 
o 

we present an EW for mixed state (1251) in 7-L 2 <8> H3 Hilbert space as follows 

4 9 

W 203 = h ® 4 + E E A A ® Qj , (49) 
i=l i=l 

where I 2 and ^3 are 2x2 and 3x3 identity matrix, respectively. A%jS are some param- 
eters, and OiS which have been defined in the previous section (see 132 al until 132 d|) and 
QjS are Hermitian operators from second party Hilbert space as following 



Qi = 


|1><1|, 


Q2 = 


|0)(0|, 


Qa 






Qa = 




+ q 2 ) 


, Q 5 = 


1 

71' 




+ q 3 ), 


Qi = 


7=2^ 


-Q2) 


, Qs = 






-Qs), 



(50 a) 

6 = ^(Q2 + Q 3 ), (506) 
9 = ^(Q2-Q 3 )- (50c) 
So, after some mathematical manipulations we get 

An = A 23 = - 2 ^- ( ^ ^-^ T ^ ) , (Bla) 

6^/2 + (x + y)2 

Al2 = A 22 = — 1 + X + y , (516) 
6 + 3(x + y) 2 
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-2v / 3-(v / 3(x + Z/) 2 + 3 v /2 + (x + ?/ )2) 

AL3 — A2I — / = , (51c) 



A 35 = A i8 = -l. (51 d) 

Obviously, we can see that the A^A = AA^ < 1 . So we have an EW candidate as 

>V 2 ®3 = h ® ia + ® Qi + C 2 ® Qs) + A 2 (Ci <g> Q 2 + 2 ® Q 2 ) + 

A 13 (Ox ®Q 3 + <D 2 ® Q x ) - (0 3 ® Q 5 + C 4 ® Qs). (52) 

Its expectation values in separable states are all nonnegative while its expectation value 
in the state (BSD rea ds 

Tr(W 203 • Pi) = ^{6 - Q^J(x-y) 2 - 3yJ(x + y)> - ^6 + 3(x + y) 2 }. (53) 

We want to show that for some values of the parameters x and y, the amount of en- 
tanglement of mixed state (125]) changes . To this aim, we have constructed a non-linear 
EW (152|) . For convenience, we consider point x = 0,y = 1 which is entangled , i.e., 
Tr(W 2( g)3 ' p) — — 1- So, in another inertial frameS" that moves with velocity v with 
respect to rest frame with 9 — ?, we have 

n = -(1 -x-y), (54a) 
r 2( 3) = ^(2 - 2x + Ay ± 3^3x2 + y 2 ), (545) 



73(4) = ^(2 + 4x - 2y ± 3^ + 3?/ 2 ), (54c) 

where Tj's are eigenvalues of the partial transpose = C/(A) p% U(AY . So, the non-linear 
EW, W 2( g,3, can be constructed in a similar way for relativistic density matrix. Its ex- 
pectation value with respect to relativistic density matrix p± is given by 



where 



Tr(HW^) = l~{y/(x - y) 2 +^(x + y) 2 +^B- 2^A+ + 2^1}, (55) 

B = 4 + llx 2 - 8xy + lly 2 , (56a) 

A = 4 + 4x 2 (l + x)(5x - 1) - Axy(Ax + llx 2 - 9) + 

y 2 (97x 2 - 16a; - 4) + 4(4 - llx)y 3 + 20y 4 . (566) 

After some calculations one arrives at 

Tr[W 2 ®z- P i} <Tr[W 2 ®s-p$}. (57) 

For example, for x = 0,y = 1 we obtain Tr(W • p±) = — 7p, this result shows that 
the measure of entanglement single spin one particle decreases when the velocity of the 
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observer increases. In Fig. 1, we have shown that the effect of the Lorentz transformation 
is to increase the region of separable states, (see Figure. 1) 

Moreover, using the geometric measure of entanglement for quantum systems defined 
in [2TJ 122], we consider entangled state ( 125]) with x = and y = 1 from the nearest 
separable state in the case x — |,y = |, which have been obtained using the convex 
optimization. We have 

, l9 cos 2 9 , , 

\\P-Ps\\ 2 = ^-, (58) 

and 

-1 



Tr[W ■ p] = — cos 0^6 -2 cos 20, (59) 
which is separable in case of 9 — 5. 



4. SUMMARY AND CONCLUSION 



In this paper, we have considered spin-momentum correlation of single spin half and 
one relativistic particle quantum states by using the NLEW, which in case of spin 
half coincides exactly with the concurrence. Then, we have shown that for single- 
particle quantum states which have been constructed based on two types of degrees 
of freedom spin and momentum, both the amount and the region of entanglement 
between spin and momentum decrease under Lorentz transformation, with respect to 
the increasing of observer velocity. Likewise, we have obtained the nearest distance 
bound of separable states from entangled density matrix and shown that it leads to 
zero in the ultrarelativistic limit. Finally, a natural question arises as how the previous 
results would generalize to the case of other mixed states. So, the calculations in this 
study are intended as a point of reference for the development of an understanding of 
the measure of entanglement in critical quantum systems based on NLEW. 



APPENDIX A 



Wigner representation for spin-| 

In Ref. [25] , it is shown that the effect of an arbitrary Lorentz transformation A 
unitarily implemented as U (A) on single-particle states is 



U(A)(\p) ® \a)) = ^^£D^WA,P))(|AP> ® 1*0), 

where 

W(A,p) = L- 1 (Ap)AL(p), 

is the Wigner rotation [8]. We can view this Wigner rotation as follows: we perform 
the Lorentz transformation L(p) on the rest frame to obtain a moving frame 1, followed 
by a transformation from frame 1 to frame 2 with A. Then we return to the rest frame 
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by further performing L _1 (Ap). This rotation of the local frame of rest is the kinematic 
effect that causes the Thomas precession. We will consider two reference frames in this 
work: one is the rest frame S and the other is the moving frame S' in which a particle 
whose four-momentum p in S is seen as boosted with the velocity v. By setting the boost 
and particle moving directions in the rest frame to be v with e as the normal vector in 
the boost direction and p, respectively, and h = e x p, the Wigner representation for 
spin- 1 particle is found as [7], 




where 




sin — n 
2 




and 



cosh a = 7 



1 



cosh 5 = — 
m 
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Figure 1. 
Figure Caption 

• Fig. 1: (Color online) Here, the space of the density matrices under Lorentz 
transformation is shown. The solid line (black triangle) shows the positivity 
condition. The hatched area with the dashed lines is the boundary of separability 
in rest frame(i.e., 9 = 0). The dotted black curve represents the boundary of 
separability density matrix under Lorentz transformation when 9 — |. For 9 = | 
leads to black triangle and positivity condition in rest frame. 
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Figure 2. 
Figure Caption 

• Fig. 2: The Tr\V\?2®2 * Pi] versus Wigner angle Vt p . Here pi = |, p 2 = 0, p% = |, 
and p 4 = h- 



